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Abstract. Lectures deal with the theory of electronic transport, in particular with the electrical 
conductivity, in systems dominated by strong electron-electron repulsion. The concept of charge 
stiffness is introduced to distinguish conductors and insulators at T = 0, but as well usual resistors, 
possible ideal conductors and ideal insulators at finite temperature. It is shown that the latter 
singular transport appears in many integrable systems of interacting fermions, the evidence coming 
from the relation with level dynamics, from the existence of conserved quantities as well as from 
numerical studies and exact results. Then, exact duagonalization approaches for the calculation of 
static and dynamical quantities in small correlated systems are described, with the emphasis on the 
finite-temperature Lanczos method applicable to transport quantities. Finally, anomalous dynamical 
conductivity within the planar single-band model is discussed in relation with experiments on 
cuprates. 

INTRODUCTION 

In the theory of strongly correlated electrons there are at present still numerous open 
problems, both regarding the basic understanding of phenomena and even more of the 
methods of calculation and sensible approximations. Among such theoretical challenges 
are also electronic transport quantities, such as electrical conductivity, the electronic 
heat conductivity, the thermopower and the Hall constant. In usual metals and semi- 
conductors the transport theory, based on the Boltzmann semiclassical approach as well 
as the quantum calculation of relevant scattering mechanisms, involving impurities and 
electron-phonon coupling, date back to the beginning of then solid state theory. The 
extention of these methods to strong disorder and strong electron-phonon coupling, in- 
troducing the phenomena of localization, polaronic transport etc., using the framework 
of general linear response theory have been in the focus of theoreticians for several 
decades. Within the same treatment also the effects of the weak electron-electron cou- 
pling has been understood, as summarized within the Landau Fermi liquid concept for 
electron transport. On the other hand, the strong electron-electron Coulomb repulsion, 
being the signature of strongly correlated systems, raises new questions of appropriate 
theoretical techiques and moreover of the basic understanding of the scattering mecha- 
nism and its effects. 

As a prominent example of the transport quantity we will mainly focus on the dynam- 
ical (optical) electrical conductivity o(co), whereby the easiest accesible experimental 
quantity is the d.c. resistivity p = l/o(0). Within the usual Botzmann theory for a weak 
electron scattering the relaxation-time approximation represents well the low-frequancy 



behaviour, 



o(to) = g /(1 + /cot), (1) 

where the relaxation time x depends on the particular scattering mechanism and is in 
general temperature dependent. In the following we will consider only homogeneous 
systems without any disorder, so the relevant processes in the solid state are electron- 
phonon scattering and the electron-electron (Coulomb) repulsion. When the latter be- 
comes strong it is expected to dominate also the transport quantities. In this case there 
appear evidently fundamentally new questions: 

a) The repulsion can transform a metal (conductor) into an (Mott-Hubbard) insulator 
at T = which radically changes transport quantities, both in the ground state as well as 
at finite temperatures T > 0. 

b) Even in the metallic phase it is not evident which is the relevant scattering process 
determining x(T) and p(T). In the absense of disorder and neglecting the electron- 
phonon coupling the standard theory of purely electron-electron scattering would state 
that one needs Umklapp scattering processes to yield finite x. That is, the relevant 
electron Hamiltonian includes the kinetic energy, the lattice periodic potential V and 
the electron-electron interaction H[ nU 

H = H kin + V+H int . (2) 

Then, in general the electronic current density j is not conserved due to Umklapp pro- 
cesses in the electron-electron scattering (where the sum ef electron momenta is nonzero, 
equal to nonzero reciprocal wavevector), i.e. [H, j] ^ 0, leading to the current relaxation 
and dissipation. The interplay of V and H int becomes however quite involved in the case 
of a strong electron repulsion. This becomes quite clear in examples of integrable tight 
binding models of interacting fermions, which have anomalous (diverging) transport 
quantities. 

c) Experiments on cuprates, manganites and other novel materials - strange metals - 
with correlated electrons question the validity of the concept of the current relaxation 
rate 1/x. For example, in cuprates only strongly frequency (and T) dependent x((0, T) 
can account for experiments in the normal state, leading to the concept of the marginal 
Fermi liquid. 

In these lectures we will consider only the physics of strongly correlated electrons. 
The standard framework for theoretical investigations in this case are microscopic tight- 
binding models for electrons (in general with spin), including the hopping term T 
(representing the effect of V and H^ in and the electron-electron interaction Hi nt 

# = - E /,,(<•;;,<•,, -//.r.) +H int = T+H int , o) 

where c ;v (cj v ) are annihilation (creation) operators for fermions on sites i on a D- 
dimensional lattice and in most cases = t is the hopping matrix element between 
nearest neighbor (n.n.) sites only. H int is a local (onsite or including only few neighboring 
sites) interaction defined for specific cases lateron. 



LINEAR RESPONSE THEORY 



The linear response theory for the electrical a.c. conductivity a(cfl) has been formulated 
by Kubo and is standard for dissipative systems [1]. Nevetheless, in correlated systems 
at T = as well as at T > in particular cases (of integrable models) one has to take 
in general into account besides the regular part o reg ((o) also the singular contribution to 
the charge dynamical response, i.e. 

a(to) = 2nelD8((i)) +o reg ((o), (4) 

where D represents the charge stiffness. The analysis in the presence of finite D > 
as well as the derivation of the proper optical sum rule within the models of correlated 
elctrons (3) requires more care and goes beyond the standard linear response formula- 
tions, hence we present it shortly below. 

We follow the approach by Kohn [2] introducing a (fictitious) flux (|) through a torus 
representing the lattice with periodic boundary conditions (p.b.c). Such a flux induces a 
vector potential A, being equal on all lattice sites. In lattice (tight-binding) models with 
a discrete basis for electron wavefunctions the vector potential A can be introduced via 
a gauge transformation (Peierls construction) 

c )s -> c )s ex V ( - ' e oA -Rj), (5) 

which effectively modifies hopping matrix elements Uj. Taking A as small and assuming 
that magnetic field does not modify within the lowest order the (local) interaction term 
Hi nt one can express the modified tight-binding Hamiltonian (3) 

H(A) = -Y,tije- ie ^c) s c is + H int * 

e 1 

« H(0) + e Q A-j + ^-A-xA = H(0)+H', (6) 

where R ;; = Rj — R, and j and x are the current and the stress tensor operators, respec- 
tively, given by 

j = i^tijllijCjgCis, 

i,j,s 

x = JT/z/R// Ri/ h c is . (7) 

Note that in usual n.n. tight-binding models x is directly related to the kinetic energy 
operator, x aa = (H kin ) a a- 

The electrical current j e is from the equation (6) expressed as a sum of the particle- 
current and the diamagnetic contribution, 

j e = -dH/dA = -e j-e 2 xA. (8) 

The above analysis applies also to an oscillating A(/) = A((o)exp(— i(ti + t) with (0 + = 
(O + z'8 with 8^0. This induces an electric field in the system E(/) = —dA(t)/dt. 



We are interested in the response of (j e )(to). Evaluating (j) as a linear response [1] 
to the perturbation H' in Eq.(6), and with A(a)) = E((o)/?(0 + , we arrive at the optical 
conductivity, 

POO 

X(to) = // ^ to+f ([j(0,j(0)]). (9) 

Complex a(cfl) = c((o) + /c"((fl) satisfies the Kramers-Kronig relation. Since %((0 — > 
00) — > 0, we get from the equation (9) a condition for &"(u> — > 00) ; 

o((o)J(o = -^(x), (10) 

which corresponds to the optical sum rule. It reduces to the well known sum rule 
for continuum electronic systems, as well as for n.n. hopping models where (x a «) = 
— {H kin )/d [3]. We can now make contact with the definition (4). From the expression 
(9) it follows 

Axcx = i lim n <a(») = ^[(W -JC'aa(O)], (11) 

and 

e 2 1 - e^P" 
<W<°) = J^ptaa(®)=4 C aa (t0), 

1 f 00 

Caa(co) = -Re/ dte m (j a (t)j a (0)). (12) 

In any finite system one can write C aa ((o) and D aoc in terms of exact eigenstates \n) 
with corresponding energies £„, 

C«a((0) = ^ £ p„|(m|j a |n)| 2 8((0-£ m +£ n ), 

Axa(co) = i[(Taa)-£pn K T l7 y 2 ], (") 

where p„ = e~P e "/Z is the Boltzmann weight of each level and Z = £„e~P e " is the 
thermodynamic sum. 



Conductors and insulators 

Ground state 

At T = the charge stiffness = D aa (T = 0) is the central quantity determining 
the charge transport. As already formulated by Kohn [2] D° can be expressed directly as 



a sensitivity of the ground state to the applied flux (|) a = eA a , 



i a 2 £o| 

2N 34>2 l<h=o- 



^d^o, 
^aa U..— n- v A ^7 



Since at T = there cannot be any dissipation and one expects o reg ((i> — > 0) = 0, we have 
to deal with two fundamentally different possibilities with respect to D° (for simplicity 
we consider an isotropic tensor D°): 

a) D° > is characteristic for a conductor or metal, 

b) D° = applies to an insulator, which can have the origin in the filled band of 
electrons (usual band insulator) or for non-filled band in the Mott-Hubbard mechanism 
due to strong electron repulsion, i.e. electron correlations. The latter situation is clearly 
of interest here. Note that the same criterion of the sensitivity to flux has been be applied 
to disordered systems, relevant to the theory of the electron localization. 

The theory of the metal-insulator transition solely due to the Coulomb repulsion 
(Mott transition) has been intensively investigated in last decades and is one of better 
understood parts of the physics of strongly correlated electrons. The emphasis has been 
on analytical and numerical studies of particular models of correlated electrons. The 
prototype and the most investigated model in this context is the Hubbard model with the 
onsite repulsion, 

Hi m = U^n^riii. (15) 

i 

Let us consider as an example the Hubbard model on a one dimensional (ID) chain 
with L sites. We first derive the result for free fermions at U = 0. The eigenvalues in the 
presence of a flux § become £^ = —2tcos(k + §) giving: 

E = -4t cos (A: + (()), k = In integer/L. (16) 

\k\<k F 

Following Eq.(14) we get 

D° = -4(0|r|0> = -sink F = -sin^, (17) 

L 71 71 2 

where n is the density of fermions (n = 2kp/n). We notice that D vanishes for an empty 
band n = 0, i.e. D° ~ tn for n — > 0, and for a filled band n = 2 where D° ~ t{2 — n) for 
n — > 2. D° is maximum at half filling where D°(n= 1) =2t/%. 

Another simple limit is U = °°. Since in this case the double occupation of sites is 
forbidden due to Eq.(15), i.e. n^n^ = 0, fermions behave effectively as spinless and the 
result is 

D°= -|sin(7M)|, (18) 

71 

and in contrast to noninteracting case D° vanishes at half filling, i.e. D° ~ t\l — n\ for 



Analytical results in ID indicate that D u = at half filling persists in the Hubbard 
model at all U > 0, whereby the dependence D(n) is between both limits U = and 
U = oo. The insulating state at half filling is a generic feature also for a wider class of 
other ID models, as the t-V model (equivalent to the anisotropic Heisenberg model), the 
t-J model etc. 

The metal - insulator transition in higher dimensions is more difficult and open subject 
and there are no exact analytical results [4]. Still it is established that in the 2D Hubbard 
model there is an insulator at half filling for arbitrary U due to the nesting of the Fermi 
surface. The insulating character is explicitly evident also within the 2D t-J model, 
which represents the large U limit of the Hubbard model. The metal - insulator transition 
has been also extensively studied in the limit of infinite dimensions [5]. One would like 
to know the behavior of weakly doped Mott-Hubbard insulator at n < 1. Numerical 
results for D on Hubbard model as well as the t-J model in 2D seem to indicate linear 
variation D« (1 — n) a with a = 1 [6], but there are also arguments for a > 1 [4]. 



Finite Temperature 

At T > there are in principle more possibilities, since we have to deal with do = 
o reg ((D -> 0) > as well as with D(T) > [7]: 

a) Normal conductor or resistor has a nonsingular response, i.e. D = 0, and a finite 
conductivity (resistivity) do > 0. 

b) D(T) > would be a signature of an ideal conductor, which in analogy to T = 
does not show any Ohm's energy dissipation, but rather a reactive response (accelera- 
tion) of the charge in an applied external electric field. In this case do is irrelevant since 
a singular D > term dominates at low (0. We will see that such a case is realized in a 
number of nontrivial integrable many-fermion models. 

c) The situation with D(T) =0 and do = would correspond to an ideal insulator, 
since there will not be no conduction in spite of T > 0. This is clearly a very unusual 
possibility, for which so far at least indications exist in integrable systems [7]. 

TRANSPORT AND INTEGRABILITY 

It is plausible that integrable systems can have anomalous transport properties. In classi- 
cal systems this observation is well established. On the other hand, the proposition that 
integrable quantum many-body systems show dissipationless transport at finite temper- 
atures is rather recent [8]. The key to this idea comes from a study of the Drude weight 
at finite temperatures D(T > 0) . Since all integrable many-body fermionic models are in 
ID, we will further restrict ourselves to ID, although some statements are more general. 



Transport and level dynamics 



It is convenient to discuss D(T) with a generalization of the original Kohn's approach 
[2] for zero temperature, by relating D(T) to the thermal average of curvatures D n 
of energy levels £„ subject to the fictitious flux (|) = eA. Since the flux induces the 
perturbation H' = we can expand the eigenenergies £«((()) up to the second order 
in 0, 

8,(40 = 8,(0) -Hn\M) + \{n\j\n) + <$> 2 £ IMMf. (19) 

Following Eq.(13)) we can then write D(T) as 

D = -£p„D„ = — £ Pn -^-|^ , (20) 

i.e. in analogy to T = case D can be expressed in terms of curvatures of levels, at T > 
weighted with their Boltzmann factors. 

Evaluating the second derivative of the free energy F as a function of flux (]), we also 
find 

In this expression the third term vanishes by symmetry and in a system as ours without 
persistent currents, d 2 F /d§ 2 has to vanish in the thermodynamic limit. Therefore in the 
L — > °o limit the charge stiffness can be alternatively represented as 

On the other hand, we know from rather recent studies that the integrability of a 
system is reflected in the distribution of level spacings P(s) [9]. For integrable sys- 
tems P(s) follows the Poisson distribution P(s) exp(-s), while in non-integrable 
systems the Gaussian orthogonal ensemble (GOE) one given by the Wigner surmise, 
P(s) oc s exp(—s 2 ). The difference is due to the existence of a macroscopic number of 
conservation laws in integrable systems that, allowing level crossing, supresses level re- 
pulsion. As for the behavior of D, we can now argue that in non-integrable systems, due 
to the level repulsion, the levels move as a function of § in an energy spacing of the or- 
der of the inverse of the many - body density of states. The latter becomes exponentially 
small with increasing L, leading to an exponential suppression of slopes de n /d§ as well 
as D(T). On the other hand, in the integrable system, the slopes due to the absence of 
level repulsion can take values of the order of one and D(T) can remain finite even in 
the thermodynamic limit. 

In non-integrable systems the assumption of the random matrix theory has also other 
consequensces. In particular, there exists a relation between off-diagonal and diagonal 



elements [10] 

(I jnm \\ n ^ m = (Unn-V(E)) 2 )=C 2 (E). (23) 

Here j nm = (n \ j \ m), v = (j nn ) and the outer brackets denote an average of matrix 
elements locally in the spectrum. This leads for a normal resistor with D(T > 0) = to 
the relation for the d.c. conductivity, 

ao = (3rce 2 £ Pn | jnm | 2 5(e n -B m ) = ^ f° ^p 2 (£)o 2 (£)ffi. (24) 

n-£m 

Oo can be thus expressed in terms of o(e), which measures the sensitivity to flux 
(velocity of level dynamics). Such an approach has been used already in the theory 
of conduction in disordered systems. 



Models 

Particle in afermionic bath 

The first model we study is a toy model of dissipation. The Hamiltonian describing a 
single (tagged) particle interacting with a bath of spinless fermions is 

H=-t h ^(e i ^ +1 d i + Hx.)-t^(cl 1 c i + H.c.) + U^d]d i c]c i , (25) 

( 

where ci{c\) are operators for spinless fermions (without charge, i.e. not sensitive to the 

flux) and di(dj) for the tagged particle on a chain with p.b.c. The interaction comes only 
through the on-site repulsion U > 0. The current j refers to the tagged particle only, 

j = -it h J £d] +l d l +H.c. (26) 

i 

and is not conserved for U > 0, i.e. [j,H] ^ 0. In this model the conductivity a(co) 
corresponds to the mobility //(to) of the tagged particle. The model (25) is integrable by 
the Bethe ansatz method in the case of equal masses, i.e. t h = t (where it is equivalent 
to the problem of a Hubbard chain in a nearly polarized state S z = S z max — 1) and non- 
integrable for unequal masses t^^t. 

In Fig. 1 we present D(T) as a function of t^/t [8], for different T and L but for fixed 
U/t = 2 at half-filled band of spinless electrons N e = L/2. Results clearly reveal: 

a) a nonmonotonic behavior of D(T) as we sweep through the integrable t n = t point 
for T > 0, 

b) a very weak dependence of D(T) on the system size L for the integrable case t n = t, 
and on contrary a rather strong one for the non-integrable ones. 

To obtain an impression of the overall behavior of D(T) as a function of T, we present 
in Fig. 2 D(T)/Dq for t^/t = 1.0, 0.5 and different size systems. 
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FIGURE 1. D(T)/Dq as a function of th/t for U/t — 2, different temperatures T and system sizes L: 
(- • -)L = 6, (- • • -)L=10, (- • • • -)L=14. 




Above results are clearly consistent with the conjecture [8] that at the integrable 
point th = t we are dealing with ideal conductance (tagged particle mobility) with 
D(r > 0) > 0, while for t h ^ t the transport is "normal" with D(T > 0) = and finite d.c. 
mobility /jq. A very relevant but so far open question is how does /jq behave close to the 
integrable point, where it diverges. One would speculate on the power law divergence 
Ho \t — th\^ y with y > 1, which is partly supported by numerical results [11]. 



t-V model 



In contrast to the system of a single particle in a bath the second model, i.e. the 
generalized t-V model, is a simplest model of a homogeneous system of interacting 
spinless fermions 

L L L 

H=-t^ (e%] +l a + H.c.) + V £ mn i+1 + W £ mm+2, (27) 

('= 1 ('= 1 ('= 1 

where the interaction is between fermions on the n.n. sites (V term) as well as next n.n. 
sites (W term). This Hamiltonian is integrable using the Bethe ansatz method for W = 
and non-integrable for W ^ 0. For W = and V <2t the ground state is metallic, while 
for V > 2t a charge gap opens at half filling n = 1/2 and the system is an insulator. For 
W = the t-V model is equivalent to the Heisenberg spin 1/2 chain as it can be shown 
using a Jordan-Wigner transformation. 

We studied numerically [7] various size systems with p.b.c. and N e = L/2 fermions 
(half-filled band). In order to avoid an artificial smoothing procedure of the discrete 
o((o) spectra of our finite size systems, it is convenient to present the integrated and 
normalized intensity I(a>) 

2T f® 21 D 

/(co) =£>* + —- / rfco'o(to'), D* = -—. (28) 
{-T) Jo+ {-T) 

We show in Fig. 3 I((£>) and D* for a non-integrable case with W ^ 0. Indeed, as 
expected for a generic metallic conductor (resistor), we find that D* scales to zero, 
probably exponentially with system size L. At the same time Go > as /((*)) approaches 
(0 = with a finite slope. These two results imply that the non-integrable system indeed 
behaves as a normal conductor at T > 0. 

In contrast to the non-integrable case the behavior of 1(d)) for an integrable system , 
as shown in Fig. 4, is totally different. In order to study the finite size dependence of the 
charge stiffness, we plot in the inset D* as a function of 1/L; the dashed lines indicate a 
3 rd order polynomial extrapolation based on the L = 8, 12, 16 site systems, suggested by 
the very good agreement obtained with the T = analytical result. 

We find that for L — > °o the extrapolated D* ^ 0. At the same time Go = g((0 — > 0) — > 
as /(co) seems to approach (0 = with zero slope. This behavior is reminiscent of a 
pseudogap. These two results indicate that the integrable system behaves as an ideal 
conductor at T > 0. Moreover we find that the normalized D* approaches a nontrivial 
finite value in the limit T — > °o, depending on V jt and on band filling, as both D and 
(-T) are proportional to |3 in this limit. 

Finally, as we mentioned in the introduction, the integrability is expected to have an 
effect also in a "Mott-Hubbard" insulating state. In the traditional scenario, the insulator 
at zero temperature is characterized by Do = and Go = 0; at finite temperatures T > 0, 
D(T) remains zero but due to thermally activated conduction, Go > 0. As we will show 
below this scenario seems to be verified in the generic non-integrable systems but not in 
the integrable ones. 

The test case is the integrable t-V model (W = 0) at half-filling n = 1/2, but in the 
insulating regime V >2t. The ground state here is insulating characterized by Dq = and 
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FIGURE 4. Integrated conductivity /(a>) for V = 1 .5f , W = 0, T = 2t , for L = 8 - 16 (exact diagonaliza- 
tion - full fines) and L = 20,24 (Lanczos method - dotted lines). Inset shows normalized charge stiffness 
D* vs. exact diagonalization (disks), T > Lanczos method (circles), analytical result (square) and 
3 rd order polynomial extrapolation from L = 8, 12, 16 (dotted line). 



a finite charge gap. Numerical results at finite temperatures [7] D*(T > 0) = /((O = 0) 
seems also to decrease exponentially with the system size scaling to zero for L — > oo. 
This precludes a possibility for ideal conductivity at T > 0. Furthermore /(to) seems to 
approach (0 = with a zero slope, showing a depletion of weight within a low frequency 
region of order (0 < t. These are characteristics of an ideal insulator, not conducting 
even at high temperatures T 3> Ao. 

In contrast, non-integrable systems of roughly the same charge gap exhibit a qual- 
itatively different behavior, /(to) approaches (0 = with finite (although small) slope, 
consistent with a small static conductivity Go > 0, consistent with the conduction due to 
electron - hole pairs, thermally activated over the charge gap. 



Transport and conservation laws 

The integrable quantum many-body systems, as the Heisenberg spin 1/2 chain or 
the one dimensional Hubbard model, are characterized by a macroscopic number of 
conserved quantities. It is natural to think that the singular transport behavior is related 
to the existence of such conserved quantities. A set of them is represented by local 
operators Q n , commuting with each other [Q ni Q m ] = and with the Hamiltonian, 
[Q n ,H] = 0. The index n indicates that the operator Q n is of the form Q n = Ef =1 g", 
where q" are local operators involving n sites around site /, on a lattice of L sites. 

The physical content of Q n is not clear in general. Nevertheless, for different models 
discussed below, the first nontrivial quantity Q3 (£>2 is the Hamiltonian) has a simple 
physical significance. It is equal or closely related to the energy current operator j E . 
Through its coupling to the current operator j it gives an alternative explanation for the 
singular D(T > 0) > 0. 

The time decay of correlations is related to conserved quantities in a Hamiltonian 
system by an inequality introduced by Mazur[12, 13] 

^ An ^2 

(29) 



t^tJo "V <Ql> 



which is based on a general positivity of spectra which hold for stationary processes. 
Here, operators are assumed hermitian, A^ = A, and <> denotes the thermodynamic 
average. The left hand side can be writen as a sum of a time-independent constant D and 
a time-dependent (decaying) function C(t), expressed in terms of eigenstates, 

<A(t)A> = Pn\ < n\A\m > | 2 + £ p n \ < n\A\m > \ 2 e '^n-£ m )t 

= D + C(t). (30) 

For time correlations < A(t)A > with a non- singular low frequency behavior the second 
term lim^oo j Jo C{t)dt goes to zero and therefore 

D>r <A ^ >2 . (3D 
r <qi> 



We can use this inequality for the analysis of c((o) related to the current-current 
correlation < j(t)j >. Charge stiffness D, Eq.(22), in this case is directly related to 
D, and we get 

D >l E £Za>!. (32 ) 
<Ql> 

In this derivation we assume that the regular part o reg ((£>) shows a non-singular behavior 
at low frequencies so that the contribution from C{t) in (29) vanishes. The latter is 
consistent with numerical simulations showing for integrable systems a pseudogap at 
small (0 and so a vanishing regular part o reg ((£> — > 0). 

In general it is difficult to evaluate the right hand side of the inequality (29) involving 
the "overlap" < AQ n >. However, in some examples presented below, this correlation 
can easily be evaluated in the case of a grand canonical trace over states, in the thermo- 
dynamic limit and for |3 — > 0. 

Before studying concrete models we construct the energy current operator j E for tight- 
binding models with n.n. hopping only, where the Hamiltonian can be presented as a sum 

L 

tf=J> M+ i. (33) 

i=\ 

where are local energy operators acting on a bond. H is a conserved quantity, hence 
the time evolution of the local hij + i can be writen as the divergence of energy current 
operator jf: 

d J^=i[H Al+1 ] = -Uf +l -jf), (34) 
and we get for the (total) energy current 

L L 

f = Y,jf = -/J>- M ,/* M +i]- (35) 
(=1 i=\ 



t-V model: Heisenberg model 

t-V model is through the through a Jordan-Wigner transformation equaivalent to the 
anisotropic Heisenberg spin (S = 1/2) Hamiltonian, provided that J x = J y , 

H = t +JyS?S y i+1 +J z Sff +1 ] . (36) 

(=i 

Here the local energy current operators jf can be expressed as 



jf = J x J y (xzy - yzx) ,•_ i )i+ 1 + J y J z (yxz - zxy) ,•_ i 1 + J Z J X (zyx - xyz) /_ i ,/+ 1 , (37) 



with (aPY-7Pa),_ M+ i = Sf^sfsl, -Cifei- 

It is straigthforward to verify that the global energy current operator j E commutes 
with the H (36) and coincides with Q3. Consistent with the notation Q2 for the Hamilto- 
nian is the sum of local operators involving two sites ^2 = ^y+b = jf involves three 
neighboring sites etc. 

The vanishing commutator [j E \H] = implies also that the energy current time 
correlations are independent of time. The non-decaying of the energy current then leads 
evidently to a singular thermal conductivity K related to the < j E (t) j E > correlation as 
well as a singular thermopower related to < j E (t)j >. 

In the language of the fermionic t-V model j E is given by 

jf = t 2 (ic] +l a-\ +h.c.)+Vjij+i(ni-i +/1/+2- 1), (38) 

where jij+\ = t(icj +l Ci + h.c.) is the local particle current. From (32) one can get an 
explicit expression for |3 — > 0, where the leading order of the high-T expansion of {j E j) 
and (jj) yields 

P 2Vyi-i»)(2i»-l)2 
- 2 l+V 2 (2n 2 -2n+l)- K ' 

We note that outside the half-filling n^ 1/2, D > implies again an ideal conductivity 
as shown before using the argument of a relation to the level dynamics. For n = 1/2 the 
inequality is however insufficient for proving that D is nonzero. Due to the electron-hole 
symmetry, this remains true even if we consider all higher order conserved quantities 
Q n . This evidently presents an open problem, since the inequality seems to be exhausted 
neglecting possible nonlocal conserved quantities. 



Hubbard model 
As above one can define a local energy operator for the Hubbard model, 
hi,i+i = -tY,(clci+L S + h.c.) + 

s 

+ \ [( w «t - \) ( n n ~\) + K-+it - \){"i+n -\)], (40) 

and 

jf = Y, t2 ( ic l+i,s c i-i,- + h.c.) - yi./', 1.,., • ./',.,. i.v )(»/. , - (41) 

where jij+i jS = (— /cj +1 s ci s + h.c.) is the particle current. By comparing the expres- 
sion for j E to Q3, we find that they coincide when the factor U/2 in (41) is replaced by 
U. One can verify that the energy current j E has a finite overlap < j E Q^ >, with the 



conserved quantity Q3. We therefore conclude that although the energy current corre- 
lations are time dependent still they do not decay to zero so that the thermal transport 
coefficients are singular. 

Using again (32) we find for the charge stiffness D of the Hubbard model, again at 

D > fifjl L,n,(l-n,)(2n_,-l) 

"2 2Z s n s (l-n s )[l + (U/2) 2 (2n 2 _ s -2n_ s + l)Y 

where n s are densities of s =f, J, fermions. Again, we get D(T) > outside half filling 
n s 7^ 1/2. On the other hand D(T) = for n s = 1/2. In the latter case the validity or 
deficiency of the result is not so clear, since the insulating D = at half filling is not 
unexpected and there are so far also no reliable analytical or numerical result which 
could be confronted with the one following from the Mazur inequality. 



NUMERICAL METHODS 

Let us turn to the discussion of higher-dimensional (as well as nonintegrable ID) corre- 
lated systems which behave as normal metals with finite d.c. transport quantities, e.g. the 
resistivity p(T). It is evident that the evaluation of Oo(T) via the linear response theory 
requires the method at T > 0. 

The absense of well-controlled analytical approaches to models of strongly correlated 
electrons has stimulated the development of computational methods. Conceptually the 
simplest is the exact diagonalization method of small systems [6]. In models of corre- 
lated electrons, however, one is dealing with the number of quantum basis states N st 
which grows exponentially with the size of the system. Within the Hubbard model there 
are 4 basis states for each lattice site, therefore the number of basis states in the Af-site 
system is N st A N . In the t-J model N st still grows as oc 3^ and within the t-V (or 
Heisenberg) model as 2 N . In the exact diagonalization of such systems one is there- 
fore representing operators with matrices N st x N st , which become large already for very 
modest values of N. 

One straightforward approach to dynamics at T > is to perform full diagonalization 
of small system, calculating all eigenstates and eigenfunctions e n , \n), required to eval- 
uate matrix elements, e.g. (n\j\m). At the present status of computers we are restricted 
to diagonalization of matrices with N st < 3000, so that reachable systems are N < 16, 
N < 12, Af < 8 for the t-V model, the t-J model and the Hubbard model, respectively. 



Lanczos method 

The helpful circumstance is that for most interesting operators and lattice models 
only a small proportion of matrix elements is nonzero within the local basis. Then, the 
operators can be represented by sparse matrices with Nst rows and at most f(N) nonzero 
elements in each row. In this way memory requirements are relaxed and matrices up to 
N s t ~ 10 8 are considered in recent applications. Finding eigenvalues and eigenvectors 



of such large matrices is not possible with standard algorithms performing the full 
diagonalization. One must instead resort to power algorithms, among which the Lanczos 
algorithm is one of the most widely known [14]. It is mainly used for the calculation of 
the ground state of a finite system, i.e. the ground-state energy Eq and wavefunction 

The Lanczos algorithm [14] starts with a normalized random vector \§q) in the 
vector space in which the Hamiltonian operator H is defined. H is applied to \§q) and 
the resulting vector is split in components parallel to |(])o), and orthogonal to it, 
respectively, 

tfN>o>=aoN>o>+fciN>i>. (43) 

Since H is Hermitian, ao = (§o\H\§q) is real, while the phase of |<f>i ) can be chosen so 
that b\ is also real. In the next step H is applied to 

H\^ l )=b\\^)+a l \^)+b 2 \^ 2 ), (44) 

where \§2) is orthogonal to and It follows also b\ = (§o\H\§i) = b\. Proceeding 
with the iteration one gets in i steps 

H\$i) = bilh-i) + ai \h) + b i+ i\h + i), 1 < i < M. (45) 

By stopping the iteration at i = M and putting the last coefficient bjy+i = 0, the 
Hamiltonian can be represented in the basis of orthogonal Lanczos functions |<|>j) as the 
tridiagonal matrix H M with diagonal elements a ; with i = 0...M, and off diagonal ones b[ 
with i = I...M. Such a matrix is easily diagonalized using standard numerical routines 
to obtain approximate eigenvalues £y and the corresponding orthonormal eigenvectors 

Ivy). 

M 

\Vj) = Y,Vji\<bi), 7 = 0.. .M. (46) 

It is important to realize that \x\fj) are (in general) not exact eigenfunctions of H, but 
show a remainder 

H\\fj) - £j\\\fj) = b M+ iv jM \§M+i)- (47) 

On the other hand it is evident from the diagonalization of Hm that matrix elements 

(yf i \H\\fj) = Zj§ ij , i,j = 0...M (48) 

are exact, without restriction to the subspace Lm- 

The identity (48) already shows the usefulness of the Lanczos method for the calcu- 
lation of particular matrix elements. As an aid in a further discussion of the Lanczos 
method we consider the calculation of a matrix element 



W k i = (n\H k BH l A\n), 



(49) 



where \n) is an arbitrary normalized vector, and A, B are general operators. One can 
calculate this expression exactly by performing two Lanczos procedures with M — 
max(&,/) steps. The first one, starting with the vector \§o) = \n), produces the subspace 
Lm = 7 = 0.. .M} along with approximate eigenvectors \\\fj) and eigenvalues £,-. 

The second Lanczos procedure is started with the normalized vector 

|$o) =A|^ )/y / (^o|AtA|^o), (50) 

and results in the subspace L M = j = 0. . .M} with approximate |\j/y) and £y. We 

can now define projectors 

fl»=I>m ^™=Ei$i)<$ii, (5i) 

i=0 ;=0 

which for m = M can also be expressed as 

M M 

ftf = ElVi)(v/l, = £|%)(%I • (52) 

By taking into account definitions (51), (52) we show that 

HP m = P m+l HP m = P M HP m , m<M. (53) 

Since in addition \n) = |<|>o) = Pol^o) and A|n) oc |<j> ) = Pol^o). by successive use of the 
first equality in (53) we arrive at 

W U = (<$>o\PmHPmH...HP m BP m H...P m HP m A\$ ). (54) 

We note that the necessary condition for the equation (54) is M > k, I. We finally expand 
the projectors according to expressions (52) and take into account the orthonormality 
relation (48) for matrix elements, and get 

M MM M 

w k i = £ ... £ £ ... £ {to\Vh){Vi \H\y il )...{Vi k _ 1 \H\v ik ) 

'0=0 i k =0j =0 j,=0 

x (¥4 \ B \Vh)(Vj, \ H \Vj,-i) ■ ■ ■ (¥71 \H\Vh) (¥jo l A l<M = 

M M 

= EE^IV/XViWyX^ylAltfroX^W- (55) 

i=0j=0 

We have thus expressed the desired quantity in terms of the Lanczos (approximate) 
eigenvectors and eigenvalues alone. 

Within the Lanczos algorithm the extreme (smallest and largest) eigenvalues £,-, 
along with their corresponding |\|/ ; ), are rapidly converging to exact eigenvalues E t 
and eigenvectors It is quite characteristic that usually (for nondegenerate states) 
M = 30 — 60 <C N st is sufficient to achieve the convergence to the machine precision of 



the ground state energy Eq and the wavefunction ^Vq), from which various static and 
dynamical correlation functions at T = can be evaluated. 

After l^o) is obtained, the g.s. dynamic correlation functions can be calculated within 
the same framework [15]. Let us consider the autocorrelation function 

C(t) = -i(^o\A\t)A\^ ) = -/(¥o|aV'( £ °-^A|¥ > (56) 
with the transform, 

C(to)=/ dte™ +t C{t) = (y |At - A\V ). (57) 

Jo (>) + +Eo-H 

To calculate C((o), one has to run the second Lanczos procedure starting with the normal- 
ized function |(j>o), equation (50). The matrix for H in the new basis L M , with elements 
(<j>;|/f |<j>y) = [HAf\ij, is again a tridiagonal one with d\ and b{ elements, respectively. Ter- 
minating the Lanczos procedure at given M, one can evaluate the C(cfl) as a resolvent of 
the Hm matrix which can be expressed in the continued-fraction form [15], 



(^o|AU|^o>_ 
(£>++E -a ! 



C(Q)) = ^ ^2 , (58) 



(Q + +E -ai 



co+ +E -d2 - . . . 

terminating with &m+i = 0, although other termination functions have also been em- 
ployed. 

The spectral function C((o) = — (l/7t)ImC((o) is characterized by frequency moments, 

/oo 
w l C(w)dw= (y \A\H-E ) l A\y ), (59) 

which are particular cases of the expression (49) for B = A', k = 0, and \n) = ^Vq). Using 
the equation (55) we can express m for / < M in terms of Lanczos quantities alone 

M 

IH = ^(^m^Mmiej-Eo) 1 . (60) 
Hence moments /j/<m are exact for given |*Pq)- 



Finite temperature Lanczos method 

The novel method for T > [16, 17] is based on the application of the Lanczos 
iteration, reproducing correctly high-T and large-(0 series. The method is then combined 
with the reduction of the full thermodynamic trace to the random sampling. 

We first consider the expectation value of the operator A in the canonical ensemble 

(A) = ^(n\e-?> H A\n) I X>|<T P ». (61) 

n=l ' n=l 



A straightforward calculation of (A) requires the knowledge of all eigenstates |^* n ) and 
corresponding energies E n , obtained by the full diagonalization of H, computationally 
accessible only for N st < 3000. Instead let us perform the high-temperature expansion 
of the exponential exp(— (3//), 

N st <*> (_R\k 

n=lk=0 K - 
N st °° (—R\k 

z = EE4r(«rt- (62) 

n=lk=0 K - 

Terms in the expansion (n \H k A\n) can be calculated exactly using the Lanczos procedure 
with M >k steps and with |<|>q) = \n) as a starting function, since this is a special case of 
the expression (49). Using the relation (55) with / = and B = 1, we get 

M 

(n\H k A\n) = £<nM) W|A|n)(e?)*. (63) 

Working in a restricted basis k < M, we can insert the expression (63) into sums (62), 
extending them to k > M. The final result can be expressed as 

N^ M 

< A > « Z- l ££e-W(n\y»)W\A\n}, 

N st M 

Z « (64) 

and the error of the approximation is of the order of |3 M+1 . 

Evidently, within a finite system the expression (64), expanded as a series in p\ 
reproduces exactly the high-T expansion series to the order M. In addition, in contrast 
to the usual series expansion, it becomes (remains) exact also for T — > 0. Let us assume 
for simplicity that the ground state |*Pq) is nondegenerate. For initial states \n) not 
orthogonal to |*Fo)> already at modest M < 60 the lowest function |\|/q) converges to 
l^o). Eq.(64) for (3 — > °° then gives, 

(A) = Y,(n\W )(W \A\n) / £ H^oXW = 

= (YolATO/OFolYo), (65) 

where we have taken into account the completeness of the set \n). Obtained result is just 
the usual ground state expectation value of an operator. 

In order to calculate dynamical quantities, the high-T expansion must be supple- 
mented by the high-frequency (short-time) expansion. The goal is to calculate the dy- 
namical correlation function at T > 0, 



(B(t)A) = Tr 



(66) 



Expressing the trace explicitly and expanding the exponentials, we get 

(B(t)A) =Z- J £ £ £ (-P + ^C-^ {n \H k BH l A\n). (67) 

Expansion coefficients in equation (67) can be again obtained via the Lanczos method. 
Performing two Lanczos iterations with M steps, started with normalized |(|)q) = \n) and 
|((>o) oc A I n), respectively, we calculate coefficients Wu following the equation (55), while 
Z is approximated by the static expression (64). Extending and resumming series in k 
and / into exponentials, we get 

N st M M 

(B(t)A)^Z- 1 £ £ £ e -P £ "/( e "- e ")(n|<)«|5|^)(^|A|n). (68) 

The computation of static quantities (64) and dynamical ones (68) still involves the 
summation over the complete set of N st states \n), which is not feasible in practice. 
To obtain a useful method, one further approximation must be made which replaces 
the full summation by a partial one over a much smaller set of random states. Such 
an approximation, analogous to Monte Carlo methods, is of course hard to justify 
rigorously, nevertheless we can estimate the errors involved. 

We consider the expectation value (A) at T > 0, as defined by the expression (61). 
Instead of the whole sum in equation (61) we first evaluate only one element with respect 
to a random state |r), which is a linear combination of basis states 

N st 

k>=£P»», (69) 

i.e. |3 m are assumed to be distributed randomly. Let us discuss then the random quantity 

A r = (r\e-^ H A\r)/(r\e-^ H \r) (70) 

We express | $ rn | 2 = l/N st + 8 rn and assume that random deviations 8„ 3 are not correlated 
with matrix elements (n\e~$ H \n) = Z n and (n\e~$ H A\n) = Z n A n . Then it is easy to show 
that A r is close to (A), and the statistical deviation is related to the effective number of 
terms Z in the thermodynamic sum, i.e. 

A r =(A) + 0(l/VZ), Z = e P £ "£z„~Tr e -P( // - £ °). (71) 

n 

Note that for T — > oo we have Z — > N st and therefore at large N st a close estimate of 
the average (71) can be obtained from a single random state. On the other hand, at finite 
T < oo the statistical error of A r increases with decreasing Z. In the FTLM we replace the 
full summation in the expression (61) with a restricted one over several random vectors 
\r),r= l,R and the statistical error is even reduced, 



A = {A) + 0(1/VRZ). 



(72) 
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FIGURE 5. The variation of Tf s with doping Cf, in the t-J on N = 18 sites with J/t — 0.3. Curves 
correspond to different choices of Z* = Z(Tf s ). 



We introduce and justify the finite-temperature Lanczos method as a method to cal- 
culate T > properties on small systems, but as well as a controlled approach to results 
at T — * in the thermodynamic limit. We also argue that choosing appropriate M, and 
the sampling R one can reproduce exact results to prescribed precision on a given sys- 
tem. The well known deficiency of the exact diagonalization methods is however the 
smallness of available lattices, hence it is important to understand the finite size effects 
and their role at T > 0. We claim that generally T > reduces finite size effects. This 
is related to the fact that at T = both static and dynamical quantities are calculated 
only from one wavefunction |*Po), which can be quite dependent on the size and on the 
shape of the system. T > introduces the thermodynamic averaging over a larger num- 
ber of eigenstates. This reduces directly finite- size effects for static quantities, whereas 
for dynamical quantities spectra become denser. 

The effect of T > can be expressed also in another way. There are several char- 
acteristic length scales in the system of correlated electrons, e.g. the antiferromagnetic 
correlation length the transport mean free path l s , etc. These lengths decrease with 
increasing T and results for related quantities have a macroscopic relevance provided 
that the lengths become shorter than the system size, e.g. l s < L where L is the linear 
size of the system. This happens for particular T > T s , where clearly T s depends also on 
the quantity considered. 

As a criterion for finite size effects we use the characteristic finite-size temperature 
Tf s . It is chosen so that in a given system the thermodynamic sum in Eq.(71) is apprecia- 
ble, i.e. Z(Tf s ) = Z* ^> 1. The finite-temperature Lanczos method is thus best suited just 
for quantum many-body systems with a large degeneracy of states, i.e. large Z at low T. 
This is the case with doped antiferromagnet and the t-J model in the strong correlation 



regime J < t. To be concrete we present in Fig. 5 the variation of Tf s with the doping 
Ch = Nh/N, as calculated from the system of N = 18 sites and J/t = 0.3. It is indica- 
tive that Tf s reaches the minimum for intermediate (optimum) doping Ch = c* h ~ 0.15. 
Away from such optimum case T/ s is larger. We claim that small T/ s and related large 
degeneracy of low-lying states are the essential features of strongly correlated system 
in their most challenging regime (so called underdoped and optimally doped regime of 
cuprates), being a sign of a novel quantum frustration. This is consistent with the exper- 
imental results on cuprates, where the optimal doping with respect to T c coincides with 
largest degeneracy, i.e. electronic entropy density s. 



T-J MODEL VS. CUPRATES 

The main motivation to understand transport in strongly correlated systems comes from 
the well established anomalous properties of high-r c cuprates, whereby in our context 
we consider only the "normal" metallic state. The prominent example is nearly linear 
in-plane resistivity p T in the normal state [4]. It is however an experimental fact 
that such a behaviour is restricted to the optimum doping regime, while deviations from 
linearity appear both in the underdoped and overdoped regimes, being still universal for 
a number of materials with a similar doping. The d.c. resistivity p is intimately related 
to the optical conductivity a (to), which has been also extensively studied and shows in 
the normal state the unusual non-Drude behaviour fitted with the marginal Fermi-liquid 
[18] behaviour x ~ 2tiX((0 + ^T). 

There is a lot of evidence that the anomalous electrical transport in cuprates is mainly 
due to strong correlations. Moreover, it seems that at least within the normal metallic 
state the physics of cuprates can be well captured within the single-band t-J model, 
the simplest model which contains the interplay of the charge propagation and spin 
(antiferromagnetic) exchange, 

H = ~t £ (c) s c is + H.c.) +/£(S ; • Sj - l -n inj ). (73) 

(«7> («7> 

Here S ; are the local spin operators interacting with the exchange parameter J. Due to 
the strong on-site repulsion the double occupion of sites is explicitly forbidden and we 
are dealing with projected fermion operators c ;v = c ;v (l —rii- s ). 

Due to small number of local degrees of freedom (/V ; = 3) the t-J model is particularly 
convenient for studies with the exact diagonalization method. Ground state properties of 
the planar model on a square lattice have been intensively investigated [6], more recently 
also the T > properties using the finite-temperature Lanczos method [17]. 

Let us go straight to results at the intermediate (optimum) doping, where we can 
reach lowest Tf s . Instead of (j(cfl) directly it is more instructive to present the current 
correlation function C(cfl), equation (12). To avoid the ambiguities with an additional 
smoothing, we plot the corresponding integrated spectra 



r(0 

7 c (to)= / C{d)d(*'. 
Jo 



(74) 




FIGURE 6. Integrated current correlation spectra /c(g>) at c/, = 3/16 for different r < t . 



In Fig. 6 we present /c(co) for the intermediate (optimal) doping = 3/16, for various 
T <t. Spectra reveal several remarkable features [17]: 

a) For T < J spectra /c(co) are essentially independent of T, at least for available 

T > Tf s . 

b) Simultaneously the slope of 7c(co < 2 t) is nearly constant, i.e. we find C(co) ~ Co 
in a wide range co < It. At the same time Co is only weakly dependent on J as tested for 
J/t = 0.2 -0.6. 

c) Even for higher T > J the differences in the slope Co, as also in Ic(°°), appear as 
less essential. 

We conclude that C(co < 2 t) ~ Co implies a simple universal form, 

o(©)=C e§ • (75) 

Such a c(co) shows a nonanalytic behaviour at CO — >• 0, starting with a finite slope. This 
is already an indication that a(co) is not consistent with the usual Drude form, but rather 
with a marginal concept [18] where the only co scale is given by T . It is also remarkable 
that the form (75) trivially reproduces the linear law p oc T as well as the non-Drude 
fall-off at CO > T. It is evident that the expression (75) is universal containing the only 
parameter Co as a prefactor. 

Experimental results and theoretical considerations are often discussed in terms of 
the co-dependent relaxation time x and the effective mass m*. These can be uniquely 
introduced via the complex a(co) and the corresponding memory function M(co), 



and 



1 M"((o) 
x(ro) ~ l+M^co)/©' 
m*(co) 2c»f ^ | M'(to) . (77) 



where m t = l/lta^ is the bare band mass. Using relations (76), (77) one can formally 
rewrite a(cfl) in the familiar Drude form 

a(a>) = ? * Che2 \ - - . (78) 
v ; ot*(q))[q) + i/t(q))] 

Employing the equation (78) we evaluate both x(to) and m*(to) from known a(cfl). It 
follow directly from the form (75), that in the regime T < J and 0)<?we can describe 
the behaviour of 1 /x with a linear law 

x" 1 =2nX(& + t,T), X-0.09, ^-2.7. (79) 

This dependence falls within the general framework of the marginal Fermi liquid sce- 
nario [18]. It should be however stressed that the asymptotic form (79) does not allow for 
any free parameter, i.e. constants X and £ are universal and independent of any model pa- 
rameters, whereas within the original marginal Fermi-liquid proposal X is an adjustable 
parameter while £ = n. 
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FIGURE 7. Sheet conductivities o(co) for various T/t, in comparison with measurements in different 
cuprates. Experimental results refer to T < 200 K. 

Let us relate obtained model results to experimetns in cuprates. Since parameters 
of the model, corresponding to cuprates, are to large extent known, i.e. J/t ~ 0.3 
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FIGURE 8. Sheet resistivities p(r) for various dopings (full lines) in comparison with measurements 
in LSCO with x = 0.15 (dotted), BSCCO (dashed), and YBCO (dash-dotted). 



and t ~ 0.4eV, the direct comparison of a((0 and p(T) can be performed without any 
additional assumptions. In 2D o is naturally expressed in terms of the universal constant 
Po = h/ e o- The corresponding 3D conductivity of a stack of 2D conducting sheets with 
an average distance c is given instead by 030 = o/c. For comparison with experiments 
we reduce the 3D measured values to the 2D conductivities for three different cuprates 
at intermediate doping, i.e. La2- x Sr x Cu04 (LSCO) with x = 0.2, Bi2sr2CaCu20 8+ s 
(BSCCO) and YBa2Cu307 (YBCO). Taking into account the uncertainty in effective 
hole doping within these materials calculated spectra cr(cfl) for Ch = 3/16 are in a 
quantitative agreement with measurements, as seen in Fig. 7.1t should be noted that also 
calculated T _1 ((fl) and corresponding parameters X and as defined by the equation 
(79) are close to the experimental ones. 

In Fig. 8 we compare calculated p(T) to the measured ones. It should be pointed out 
that there is a restricted T window where a comparison could be made since Tf s ~ 450K, 
whereby at T ~ Tf s also finite-size effects start to influence our analysis. Nevertheless, 
for the intermediate doping Ch ~ 0.2 our p(T) results match quantitatively well experi- 
mental ones for cuprates with comparable hole concentrations, i.e. for BSCCO, YBCO 
and LSCO with x = 0A5. 



Origin of anomalous conductivity 

There is so far no accepted explanation for the origin of the anomalous marginal-type 
charge dynamics in cuprates. From our analysis of the t-J model at the intermediate 
doping it seems to be quite generic feature of frustrated regime of correlated electrons. 
In the following we try to argue that there is in general a relation between the marginal- 
type o(co) and the overdamped character of quasiparticle excitations, also observed 



experimentally in cuprates via the angle-resolved photoemission spectroscopy (ARPES) 
as well in model studies. 

Let us perform for C(cfl) the decoupling in terms of single-particle (electron) spectral 
functions A(k, co) neglecting possible vertex corrections, i.e. 

C ( ro ) = V^ (V * )2 / ^ / /(- ro/ )/(« / -«M(k,a) / )A(k,( / -a)), (80) 

where /(co) are Fermi functions and vjj are (unrenormalized) band- velocity components. 
While such an approach is usual in weak scattering problems, qualitative features could 
be reasonable also for strongly correlated electrons. We assume also that A(k, co) for 
quasiparticles close to the Fermi energy (co = 0) can be generally presented as 

A(k,Q)) = - 7 Zk ^, (81) 

where quasiparticle parameters Zk, I\, £k can still be dependent on co and T. In order to 
reproduce the marginal Fermi liquid form, Eqs.(78),(79), of a(cfl) we have to postulate 
an analogous form also for the quasiparticle damping, i.e. T = y( jco| +£,T) (consistent 
with recent ARPES experiments in BSCCO), but we neglect the k dependence of T and 
Z. 

For (0 < (0*, whereby the cutoff co* appears due to the effective bandwidth or some 
other characteristic quasiparticle scale, only the behavior near the Fermi surface should 
be important. Replacing in Eq.(80) the k summation with an integral over £ with a slowly 
varying density of states 3\£ (e) one can derive 

/ rfeA(e,a/)A(e,a/-a>) = % (82) 
J n co z + 1 (co, co'j^ 

where f(co, co') = r(co / ) + r(co / - co). We are thus dealing with a function C(co), 

C(co) =C I ddf(-d)f{d-u) r J}^L 2 , (83) 

J C0 Z + 1 (CO, CO J 

depending only on the ratio (o/T, and on the parameters y,^. In the theory of Fermi 
liquids co and T variation are tightly bound so ^ should also be quite universal or at least 
quite restricted in range, whereby ^ ~ 71 is usually used [18]. 

It is evident from Eq.(83) that for y< 1 we recover C(co) strongly peaked at co = 
and consequently a(co) T / (co 2 + T 2 ) with T = 2r(co/2, T). This is just a generalized 
Drude form with l/x(co) = 2r(co/2). o(co) is in this case directly determined by the 
single-particle relaxation T. 

This relation is however not valid when one approaches the regime of overdamped 
QP excitations y ~ 1 or more appropriate ~ 1. We present in Fig. 9 results for C(co) 
for several y fixing ^ = n. While for y < 0.2 still a pronounced peak shows up at co ~ 0, 
C(co) becomes for y > 0.3 nearly constant or very slowly varying in a wide range of 
co/r. For large y we find even C(0) < C(co — > °°), approaching for y 3> 1 the ratio 
C(0)/C(~) = i/t 
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FIGURE 9. Current-current correlation spectra C(co) vs. co/r for various y at fixed ^ = n. 



The main message of the above simple analysis is that for systems with overdamped 
quasiparticle excitations of the marginal Fermi-liquid form the expression Eq.(75) de- 
scribes quite well a(co) for a wide range of parameters. It should be stressed that 
nearly constant C(co < CO*) also means that the current relaxation rate 1/x* is very fast, 
1/x* ~ CO* 3> 1/x, i.e. much faster than the conductivity relaxation scale apparent from 
Eq.(79) where 1/x °c T is determined solely by thermodynamics. 



CONCLUSIONS 

The theory of electron transport in models and materials dominated by strong electron- 
electron repulsion is still in its infancy. Present lectures try to describe some of novel 
phenomena in transport, which seem to be particular to correlated systems. At the 
same time it is evident that there are even more open problems which wait for proper 
theoretical understanding and feasible methods for systematic study: 

a) The singular transport in nontrivial integrable models of interacting fermions reveal 
an evident question which are relevant scattering processes which determine transport 
quantities. Coulomb repulsion and Umklapp processes are in general clearly not enough 
to lead to current relaxation and dissipation. As shown in integrable systems transport 
quantities remain singular. It is expected that in general any deviation from integrability 
will lead to a normal resistence, however there exists no reliable approach to transport 
in the vicinity of integrable points. 

b) Even in integrable ID systems there are several open questions. Although the re- 
lation to conserved quantities is very appealing and more operational it is questionable 
whether it can explain generally the existence of ideal conductors, e.g. the ideal conduc- 
tor in the simplest cases as the t-V model at half filling, the particle in the fermionic bath 
etc. 



c) The existence of very unusual ideal insulator has not been proven nor disproven 
yet. 

d) For higher dimensional systems there is definite need for analytical methods which 
would allow the evaluation of (at least qualitatively correct) transport quantities. Re- 
cently developed numerical methods for finite correlated systems, in particular the finite- 
temperature Lanczos method, have proven to be very useful. Their weakness is clearly in 
the smallness of reachable systems. Fortunately, the physics in several relevant models 
appears to be quite local in character so that e.g. mean free path is very short even at 
moderately low temperatures. 

e) The main motivation to study correlated systems comes from the observation of 
anomalous conductivity c((fl) in high-r c cuprates as prominent example of materials 
with strongly correlated electrons. Numerical studies as well as analytical relation with 
spectral functions with overdamped quasiparticles confirm the observed behavior de- 
scribed within the marginal Fermi liquid scenario. These materials seem to follow a 
novel quantum diffusion law where relaxation is dominated by thermodynamics only. 
However its origin is not well understood. Also it is not clear how general such scenario 
might be and whether (why) it breaks down at low enough temperatures. 
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